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The motion of fluid particles as they are pushed 
along erratic trajectories by fluctuating pressure gra- 
dients is fundamental to transport and mixing in tur- 
bulence. It is essential in cloud formation and atmo- 
spheric transport|0 g|, processes in stirred chemical 
reactors and combustion systems|&, and in the indus- 
trial production of nanoparticlesQ. The perspective 
of particle trajectories has been used successfully to 
describe mixing and transport in turbulence^, §], but 
issues of fundamental importance remain unresolved. 
One such issue is the Heisenberg-Yaglom prediction 
of fluid particle accelerations ^ based on the 1941 
scaling theory of Kolmogorov[|8[ |9p (K41). Here we 
report acceleration measurements using a detector 
adapted from high-energy physics to track particles 
in a laboratory water flow at Reynolds numbers up to 
63,000. We find that universal K41 scaling of the accel- 
eration variance is attained at high Reynolds numbers. 
Our data show strong intermittency — particles are ob- 
served with accelerations of up to 1,500 times the ac- 
celeration of gravity (40 times the root mean square 
value). Finally, we find that accelerations manifest the 
anisotropy of the large scale flow at all Reynolds num- 
bers studied. 

In principle, fluid particle trajectories are easily mea- 
sured by seeding a turbulent flow with minute tracer par- 
ticles and following their motions with an imaging sys- 
tem. In practice this can be a very challenging task since 
we must fully resolve particle motions which take place 
on times scales of the order of the Kolmogorov time, 
t v = (v/e) 1 / 2 where v is the kinematic viscosity and e 
is the turbulent energy dissipation. This is exemplified in 
Fig. nl which shows a measured three-dimensional, time 



resolved trajectory of a tracer particle undergoing vio- 
lent accelerations in our turbulent water flow, for which 
t v = 0.3 ms. The particle enters the detection volume on 
the upper right, is pushed to the left by a burst of acceler- 
ation and comes nearly to a stop before being rapidly ac- 
celerated (1200 times the acceleration of gravity) upward 
in a cork-screw motion. This trajectory illustrates the dif- 
ficulty in following tracer particles — a particle's acceler- 
ation can go from zero to 30 times its rms value and back 
to zero in fractions of a millisecond and within distances 
of hundreds of micrometers. 

Conventional detector technologies are effective for 
low Reynolds number flows Jl0|, pT[ ], but do not provide 
adequate temporal resolution at high Reynolds numbers. 
However, the requirements are met by the use of silicon 
strip detectors as optical imaging elements in a particle 
tracking system. The strip detectors employed in our ex- 
periment (See Fig. [|a) were developed to measure particle 
tracks in the vertex detector of the CLEO III experiment 
operating at the Cornell Electron Positron Collider [|l2i. 
When applied to particle tracking in turbulence (See 
Fig. ||b) each detector measures a one-dimensional pro- 
jection of the image of the tracer particles. Using a data 
acquisition system designed for the turbulence experi- 
ment, several detectors can be simultaneously read out at 
up to 70,000 frames per second. 

The acceleration of a fluid particle, a + , in a turbulent 
flow is given by the Navier-Stokes equations, 

a+ = - — + iA7 2 u (1) 
P 

where p is the pressure, p is the fluid density, and u is the 
velocity field. In fully developed turbulence the viscous 



damping term is small compared to the pressure gradient 
term[pj[ p"4| ] and therefore the acceleration is closely re- 
lated to the pressure gradient. 

Our measurement of the distribution of accelerations is 
shown in Figure ||, where the probability density func- 
tion of a normalized acceleration component is plotted 
at three Reynolds numbers. All of the distributions have 
a stretched exponential shape, in which the tails extend 
much further than they would for a Gaussian distribution 
with the same variance. This indicates that accelerations 
many times the rms value are not as rare as one might ex- 
pect, i.e., the acceleration is extremely intermittent. The 
acceleration flatness, shown in the inset to Fig. 0, char- 
acterizes the intermittency of the acceleration, and would 
be 3 for a Gaussian distribution. These flatness values 
are consistent with direct numerical simulation (DNS) at 
low Reynolds number[[l4]| and exceed 60 at the highest 
Reynolds numbers. 

The prediction by Heisenberg and Yaglom for the vari- 
ance of an acceleration component based on K41 theory 
is 

(am) = a e 3/2 ^ 1/2 %, (2) 

where ao is a universal constant which is approximately 1 
in a model assuming Gaussian fluctuations [[| ^|, [l^, O], 
However, DNS has found that ao depends on e. Con- 
ventionally this is expressed in terms of the Taylor mi- 
croscale Reynolds number, R\, which is related to the 
conventional Reynolds number by R\ = (15RC) 1 / 2 and 
is proportional to e 1 / 6 . Using this notation, DNS results 
indicate ao ~ R 1 ^ 2 for R\ < 250jjl4]], with a tendency to 
level off as R\ approaches 470[]16|]. 

Our measurement of the Kolmogorov constant ao is 
shown in Fig. Q for eight orders of magnitude of scal- 
ing in acceleration variance. We find ao to be anisotropic 
and to depend significantly on the Reynolds number. The 
ao values for both components increase as a function of 
Reynolds number up to R\ » 500, above which they are 
approximately constant. The trend in ao is consistent with 
DNS results in the range 140 < R x < 4700 [T| 0. 
However, the constant value of ao at high Reynolds num- 
ber suggests that K41 scaling becomes valid at higher 
Reynolds numbers. Weak deviations from the K41 scal- 
ing such as the ao ~ i?° 135 prediction of the multifractal 
model by Borgas[|T8|| cannot be ruled out by our measure- 
ments. 



The acceleration variance is larger for the transverse 
component than for the axial component at all values of 
the Reynolds number. This is shown in the inset to Fig. Q 
where the ratio of the Kolmogorov constants for the ax- 
ial and transverse acceleration components is plotted as 
a function of Reynolds number. The anisotropy is large 
at low Reynolds number and diminishes to a small value 
at R\ = 970. This observation tends to confirm recent 
experimental results which indicate that anisotropy may 
persist to much higher Reynolds numbers than previously 

believed[|i|[!(J. 

In summary, our measurements indicate that the 
Heisenberg-Yaglom scaling of acceleration variance is 
observed for 500 < R\ < 970. At lower Reynolds 
number, our measurements are consistent with the anoma- 
lous scaling observed in DNS [[14), |l6[] . Our measurements 
show that the anisotropy of the large scales affects the 
acceleration components even at R\ w 1000. It is im- 
possible to say on the basis of these measurements if 
the anisotropy will persist as the Reynolds number ap- 
proaches infinity. We found the acceleration distribution 
to be very intermittent, with extremely large accelerations 
often arising in vortical structures such as the one shown 
in Fig. |. 

Our results have immediate application for the develop- 
ment of Lagrangian stochastic models, some of which use 
ao directly as a model constant. These models are being 
developed and used to efficiently simulate mixing, par- 
ticulate transport, and combustion in practical flows with 
varying Reynolds numbers]^ [2TJ 53], Our research also 
has surprising implications for everyday phenomena. For 
instance, a mosquito flying on a windy day (wind speed 
1 8 km/h and an altitude of 1 meter) would experience an 
rms acceleration of 15 m/s 2 . But given the extremely in- 
termittent nature of the acceleration, our mosquito could 
expect to experience accelerations of 150 m/s 2 (15 times 
the acceleration of gravity) every 15 seconds. This may 
explain why, under windy conditions, a mosquito would 
prefer to cling to a blade of grass rather than take part in 
the roller coaster ride through the Earth's turbulent bound- 
ary layer^]. 
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Figure 1: Measured Particle Trajectory The 3- 
dimensional time-resolved trajectory of a 46 micrometer 
diameter particle in a turbulent water flow at Reynolds 
number 63,000 (R\ = 970). A sphere marks the mea- 
sured position of the particle in each of 300 frames 
taken every 0.014 ms (« r^/20). The shading indicates 
the acceleration magnitude, with the maximum value of 
12,000 m/s 2 corresponding to approximately 30 standard 
deviations. The turbulence is generated between coaxial 
counter-rotating disks [|4| ^5|] in a closed flow chamber of 
volume 0.1 m 3 with rotation rates ranging from 0.15 Hz 
to 7.0 Hz, giving rms velocity fluctuation u in the range 
0.018m/s < u < 0.87m/s. Measurements are made in 
an 8 mm 3 volume at the center of the apparatus where 
the mean velocity is zero and the flow is nearly homo- 
geneous but not isotropic. As a result of a mean stretch- 
ing of the flow along the propeller axis the rms fluctu- 
ations are 1/3 larger for the transverse velocity compo- 
nents than for the axial component. The energy dissipa- 
tion was determined from measurements of the transverse 
second order structure function and the Kolmogorov re- 
lation D NN = fCx (er) 2/3 with d = 2.13[[||. The 
dissipation was found to be related to the rms velocity 
fluctuation by e = u 3 /L with an energy injection scale 
L = (71 ± 7) mm. Using the definition of the Taylor mi- 

1/2 

croscale Reynolds number R\ = {lhuL/v) 1 the range 
of Reynolds numbers accessible is 140 < R\ < 970, (in 
terms of the classical Reynolds number 1300 < Re < 
63, 000). At the highest Reynolds number the system is 
characterized by Kolmogorov distance and time scales of 
f] = 18 jtxm and t v = 0.3 ms, respectively. 
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Figure 2: APPARATUS (a) Schematic representation of the 
CLEO III strip detector[[l2|], in which grey bars indicate 
sense strips which collect charge carriers freed by optical 
radiation. The 511 strips allow measurement of the one 
dimensional projection of the light striking the detector. 
The detector may be read out 70,000 times per second, 
(b) A combination of lenses (LI, L2, L3, L3') is used 
to image the active volume onto a pair of strip detectors 
which are oriented to measure the x and y coordinates. 
Another detector assembly may be placed on the opposite 
port (LI') to measure y and z. The flow is illuminated by 
a 6 W argon ion laser beam oriented at 45° with respect to 
the two viewports. The optics image (46 ± 7) /im diam- 
eter transparent polystyrene spheres which have a density 
of 1.06 g/cm 3 . Particle positions are measured with ac- 
curacy 0. 1 strips, corresponding to 0.7 /im in the flow. 




-20 20 

a/<a 2 > 1/2 

Figure 3: ACCELERATION DISTRIBUTION. Prob- 
ability density functions of the transverse accelera- 
tion normalized by its standard deviation at different 
Reynolds numbers. The acceleration is measured from 
parabolic fits over 0.75 t, ; segments of each trajec- 
tory. The solid line is a parameterization of the high- 
est Reynolds number data using the function P(a) = 
Cexp (-a 2 / ((1 + \a/3/a\ 7 ) a 2 )), with (3 = 0.539, 7 = 
1.588, a — 0.508 and the dashed line is a Gaussian distri- 
bution with the same variance. The inset shows the flat- 
ness of the acceleration distribution, ((a 4 ) / (a 2 ) 2 , evalu- 
ated using 0.5 t v parabolic fits) as a function of R\. 
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Figure 4: a AS A FUNCTION OF R\. Open red circles 
indicate a transverse component and open red squares the 
axial component of the acceleration variance. DNS data 
is represented by blue triangles [|l4|] and green circles flip]]. 
The error bars represent random and systematic errors in 
the measurement of the acceleration variance. There is 
an additional uncertainty of 15% in the overall scaling of 
the vertical axis for the experimental data due to the un- 
certainty in the measured value of the energy dissipation. 
The degree to which the 45 /im diameter tracer particles 
follow the flow was investigated by measuring the accel- 
eration variance as a function of particle size and density. 
The results, to be published elsewhere, confirm that the 
acceleration variance of the 45^m particles is within a few 
percent of the zero particle size limit. The inset shows the 
ratio of the ao values for transverse and axial components 
of the acceleration. 
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